Numerical test of the Cardy- Jacobsen conjecture in the site-diluted Potts model in three dimensions. 
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We present a microcanonical Monte Carlo simulation of the site-diluted Potts model in three dimensions with 
eight internal states, partly carried out in the citizen supercomputer Ibercivis. Upon dilution, the pure model's 
first-order transition becomes of the second-order at a tricritical point. We compute accurately the critical 
exponents at the tricritical point. As expected from the Cardy-Jacobsen conjecture, they are compatible with 
their Random Field Ising Model counterpart. The conclusion is further reinforced by comparison with older 
data for the Potts model with four states. 
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I. INTRODUCTION 

When two ordered phases compete, even a tiny amount 
of disorder is significant. Consider, for instance, the anti- 
ferromagnetic insulator LaCuC>4. A small La «-» Sr substi- 
tution turns it into a high-temperature superconductor. Also 
for colossal magnetoresistance oxides the importance of the 
combination of phase coexistence and chemical disorder has 
been emphasized ^ 

These examples suggest a simple, yet general question: 
which are the effects of quenched disorder on systems that 
undergo a first-order phase transition in the ideal limit 
of a pure sample? (quenched disorder models impuri- 
ties that remain static over experimental time-scales 2 ^. In 
fact, this question has been relevant in a large number of 
physical contexts. A non-exhaustive list includes nanoscale 
ferroelectricity, 3 tilt ordering, 4 ferroelectric thin films,™ ran- 
dom block copolymers, 7 ferroelectric nanodisks, 8 topological 
phases in correlated electron systems, 9 effects of multiplica- 
tive noise on electronic RLC circuits 10 and surface wavesf^E] 

Unfortunately, only for two spatial dimensions (D = 2) we 
have a good understanding of the effects of quenched disor- 
der on phase-coexistence: the slightest concentration of im- 
purities switches the transition from first-order to second- 
orderPQa 

In D — 3 we lack a general description. One should con- 
sider two different possibilities: disorder may couple ei- 
ther to the order parameter, as in the Random Field Ising 
Model (RFIM)J 15 I 16 I or it may couple to the energy, as in 
the disordered Potts modelP-^ In both cases, quenched dis- 
order is unreasonably efficient at softening the transition. It 
has been surprisingly difficult to show that the transition ac- 
tually remains of the first-order for some range of impurity 
concentration£3H2ni 

Actually, the Cardy and Jacobsen conjecture relates the two 
types of disorder by means of a mapping between the RFIM 
and the disordered Potts modelPlTie conjecture reads as fol- 
lows. Consider a ferromagnetic system undergoing a first or- 
der phase transition for a pure sample. 21 Let T be the tem- 
perature while p is the concentration of magnetic sites (see 
the generic phase diagram in Fig. [TJ. A transition line, T c (p) 




FIG. 1. (Color online) Phase diagram of the three dimensional di- 
luted Potts model for Q > 3. For small dilutions we have a first order 
phase transition line which ends up in a tricritical point [at (T t ,p t )], 
and below this tricritical point, the phase transition line continues as 
a second order one. PM and FM denote a paramagnetic and ferro- 
magnetic phase, respectively. 



separates the ferromagnetic and the paramagnetic phases in 
the (T,p) plane. In D — 3 a critical concentration is expected 
to exist, 1 > p t > 0, such that the phase transition is of the 
first-order for p > p t and of the second order for p < p t (at 
pt one has a tricritical point). When p approaches p t from 
above, the latent-heat and the surface tension vanish while the 
correlation-length diverges. The corresponding critical expo- 
nents can be obtained from those of the RFIM (see below). 

However, the Cardy-Jacobsen mapping relates two prob- 
lems unsolved in D = 3. In particular, the RFIM (the suppos- 
edly well-known partner in the conjecture) suffers from severe 
inconsistencies between analytical, experimental and numeri- 
cal work. On the experimental side, mutually inconsiste nt re- 
sults for the correlation-length exponent v were obtainedjSISl 
due to the uncertainties in the parameterization of the scat- 
tering line shape. Also, the estimate of the anomalous di- 
mension tj violates hyperscaling bounds. 22 Numerical deter- 



2 



minations of exponent v are scattered on a wide rangeJ2K3 
and hyperscaling-violating results have been reported.^ The 
order-parameter's critical exponent j3 ~ 0.01 is so small (yet, 
see Ref. 23 ), that it has even been conjectured that the transi- 
tion could be of the first order.tMLZl 

On the other hand, the investigation of the disordered Potts 
model has been mostly numerical up to now. In the conven- 
tional approach, one averages over disorder the free-energy at 
fixed temperatureWil works nicely for the second-order part of 
the critical line r c (p)j23S3but the first-order piece is plagued 
by huge sample-to-sample fluctuations of the specific-heat 
or the magnetic susceptibility. 40 Fortunately, these wild fluc- 
tuations can be avoided by averaging over disorder the en- 
tropy obtained from microcanonical Monte Carlo, 44 at fixed 
energy }^ We investigated in this way the site-diluted Potts 
model with Q = 4 states. A delicate extrapolation to infinite 
system size showed that p t < 1 . Unfortunately, the relevance 
of the RFIM universality class for the tricritical point (the core 
of Cardy and Jacobsen conjecture) could not be addressed up 
to now. 

Here we show that the Cardy-Jacobsen conjecture is veri- 
fied to high numerical accuracy in the site-diluted Potts model 
with Q = A and 8 states. This results follows from a finite size 
scaling analysis of old Q— 4 data 19 and new, extensive Monte 
Carlo simulations for 2 = 8, partly carried out in the Ibercivis 
citizen supercomputer. 45 Our analysis benefits from a recent 
computation of the RFIM critical exponents, 24 that also ex- 
ploits the redefinition of the disorder average.^ 9 -' 

In Section |ll] we summarize the main implications of the 
Cardy-Jacobsen conjecture and define our specific model. Our 
methodology, including details on simulation and statistical 



analysis, is presented in Section III In Section [IV] we present 
our main numerical evidences for the validity of the conjec- 
ture. We give our conclusions in Section [V] Finally in the 
appendix we describe how the Control Variates technique im- 
proves the determination of some important quantities. 



H. THE CARDY-JACOBSEN CONJECTURE 

Specifically, we consider the D — 3 site-diluted Potts model 
with Q internal states. 17 The spins, <7/= l,...,Q, occupy the 
nodes of a cubic lattice of linear size L, with periodic bound- 
ary conditions. Each spin interacts with its nearest neighbors 
through the Hamiltonian 



(iJ) 



(1) 



The quenched randomness is represented by the occupation 
variables e, = 0, 1 (e, = 1 means that the 2-th spin is present). 
We choose each e, independently, setting £,■ = 1 with proba- 
bility p. Each specific disorder realization is called a sample. 
The pure system is recovered for p—l, where it undergoes a, 
generally regard ed as very strong, first-order phase transition 
for Q > 3 K0 | 44 | s jj 0W m Fig. Ill the full phase diagram of 
this model. 

The Cardy and Jacobsen mapping relates the large-Q limit 
of the disordered Potts model to the RFIM. 13 At the tricritical 



point p t of the Potts model, we encounter three relevant scal- 
ing fields (see, e.g., Ref. [46b . The dilution field lies along 
the critical line T c (p). We name its scaling dimension y p . 
The thermal scaling field has dimension yj, and is respon- 
sible for the ferromagnetic transition when varying the tem- 
perature. Finally, the magnetic scaling field is related to an 
external magnetic field in Eq. ([T|. The mapping to the RFIM 
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where v 1 ^ 1111 is the correlation-length exponent, 47 rj RF1M 
is the anomalous dimension, while is the hyperscaling- 
violations exponent! 1 ^ Furthermore, the exponent of the sur- 
face tension verifies a modified Widom law: jit = D — — 1 . 
Cardy and Jacobsen predicted as well that, upon approaching 
the tricritical point p,, the latent heat in the diluted Potts model 
vanishes with the same exponent p RF1M that rules the vanish- 
ing of the order parameter in the RFIM. 



III. METHODOLOGY 
A. The Microcanonical Ensemble 

For the simulation of the model described by Eq. ([TJ we 
have used an extended microcanical method which is suitable 
to study the first order part of the transition lineF 4 -' 

We will briefly review the main facts of this simulation ap- 
proach. Using a mechanical analogy, each spin is comple- 
mented with a conjugated momenta. The total energy is the 
sum of a kinetic term, -J€ (the halved sum of the squared mo- 
menta) and the potential energy, namely the spin Hamiltonian 
of Eq. Q. 

We consider the microcanonical ensemble, where the en- 
ergy (kinetic plus potential) is kept fixed to the total value Ne, 
where N = is the total number of spins. The momenta 
can be explicitly integrated out. The entropy density s{e) and 
the microcanonical weight co(e,N; {<T,}) turn out to be 



exp[Ns(e,N)} = ^,J M , £ <o(e,N; {cr ; }) ; 



NT(N/2) 



{*} 



co(e,7V;{(J,}) = 



X = Ne- Jf spin . 



(4) 

(5) 
(6) 



The role of the Heaviside step function in Eq. |5]l is preventing 
the kinetic energy from becoming negative. 

The Monte Carlo simulation of the weight in Eq. |4| is 
straightforward. Both M etropolis and cluster methods are fea- 
sible and efficient! 19 ! 44 1 In the present work we have used the 
Swendsen-Wang algorithm 44 (see Ref.[l9]for implementation 
details). One obtains in this way mean-values at fixed e that 
will be denoted ((■ ■ ■ )) e . 
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FIG. 2. (Color online) Example of Maxwell construction (data from 
a single sample of a Q = 8 Potts model in three dimensions, with 
L = 24 and p = 0.95). The horizontal line corresponds to the in- 
verse critical temperature, obtained through Maxwell's equal-area 
rule, Eq. {TT}. Consider the region limited by the horizontal line 
J3 = j3 c and the curve (J3} e . The (negatively signed) area in the re- 
gion e° < e < e* equals the absolute value of the (positively signed) 
area in the region e* < e < e d . 



Consider the equation 

P(e)=P, or (single sample) (j8/ e .}) e = /3 



(10) 



In normal situations, /3 (e) is monotonically decreasing with 
e, so that Eq. ( fTO) has a unique solution. However, at phase- 
coexistence /3 (e) is no longer monotonically decreasing, see 
Fig. [2] Therefore, Eq. ( fTO) has three important solutions, 
named e°, e*, and e d (e° < e* < e d ): 

• The rightmost root of ( fTO] ), e d ^, corresponds to the 
"disordered phase". 

• The leftmost root of ( fT0| ), e° ^ , corresponds to the "or- 
dered phase". 



• The second rightmost root of <JX0J>, e* L ^ is a saddle-point 
among the two phases. 

Note that these three solutions do depend on L, although we 
shall not explicitly indicate it unless necessary. 

We compute the inverse critical temperature j3 c from the 
equal-area rule: 



de (/3(e)- ft) 



(11) 



A particularly important mean-value comes from the 
Fluctuation-Dissipation relation 



ds 
d~e 



where 



N-2 



Ne - Jf s P m 



(7) 



(8) 



On the view of Eq. (FTl, it might be inspiring to think of (j3) e 
as the inverse-temperature corresponding to energy density e. 
The connection between the canonical and the microcanonical 
ensembles is discussed in Ref. 48 Finally, our main observ- 
able will be j3(e), defined as 



P(e) = <P)e 



(9) 



where the overline stands for the disorder-average as com- 
puted at fixed e. 



B. The Maxwell construction 

A standard way of studying phase-coexistence in a micro- 
canonical setting is the Maxwell construction. This allows to 
compute from the curve /3(e) several important magnitudes: 
the critical inverse temperature j3 c , the energies of the two co- 
existing phases and the surface tension. Furthermore, one may 
apply the very same method to the sample dependent (j5{ £ .}) e , 
as shown in Fig. [2] We follow the numerical methods de- 
scribed in Refs. 44 and [19] We briefly summarize them now, 
for the sake of completeness. 



see Fig. [2] Note that the j3 c computed from Eq. ( fTT) does 
depend on the system size. In fact, in the thermodynamic 
limit, Eq. ( fTT| is a mere consequence of the continuity of the 
free-energy density (as a function of temperature) at the phase 
transition. In fact, recall that the free energy density can be 
expressed in terms of the inverse temperature and of the in- 
ternal energy and entropy densities: / = e — s/fi . Now, if we 
recall Eq. |7]), we see that the equality of the free-energy den- 
sities of the ordered and the disordered phases at the critical 
temperature can be recast as 



Pc(e%-el)=s(e%)-s(el) 



= / "de/3(e). 



(12) 
(13) 



This textbook reasoning can be extended to the more compli- 
cated case of a finite-system. In fact, it is easy to show, see 
Refs. |44] and |49j that Eq. ( fTT) is identical to the criterion of 
equal-height in the energy histogram. 50 Such a finite-system 
indicator of the critical temperature suffers from finite-size 
corrections of order ~ 1 /L D ^ 

Once we know J3 C , we may compute the latent-heat as 



Finally, the surface-tension, L, is calculated as 

N 



L{L) = 



2L D ~ l 



(14) 



(15) 



Note that, in order to compute integrals such as the one in 
Eq. ( fTTj ), we interpolate j5 (e) (which is numerically computed 
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FIG. 3. (Color online) Maxwell construction, see Eq. {TT}, as ob- 
tained for the sample-averaged j3(e). Data for L = 48 and several 
values of the spin concentration. The transition becomes smoother 
as p decreases (from bottom to top). In fact, for p = 0.8 the Maxwell 
construction can no longer be done (because the corresponding /3 (e) 
is monotonically decreasing with e). 



over a grid in the e-line), through a cubic spline. Statistical er- 
rors are computed using a jackknife method (see e.g. Ref.l44b. 
In the case of the sample-averaged j3(e), the jackknife blocks 
are formed from the microcanonical mean-values obtained on 
the different samples. On the other hand, when one performs 
the Maxwell construction for a single sample as in Fig. [2} the 
jackknife blocks are formed from the Monte Carlo history. 

It is interesting to compare the curves /3 (e) for fixed L= 48, 
as the disorder increases (i.e. as p decreases), see Fig. [3] In the 
limit of a pure system, p — l,P(e) displays the expected cusps 
and steps for a system with well developed geometric and con- 
densation transitions. 52 As soon as the system becomes disor- 
dered, the transition becomes smoother: both the latent heat, 
see Eq. ( |14) , and the surface tension, Eq. (15) , are sizably 
smaller for p = 0.95 than for p = 1 . This trend is maintained 
for decreasing p, to the point that the phase transition is clearly 
of the second order at p — 0.8 (for that dilution, j3 (e) is mono- 
tonically decreasing with e). We note as well that the curve 
j3(e) for p < 1 is remarkably featureless, specially if com- 
pared to its p = 1 counterpart. Actually, geometric transitions 
are also found for individual samples at p = 0.95. However, 
the energies at which this singular behavior arise depend on 
the considered sample, which results in a smooth averaged 
P(e). 



C. Finite Size Scaling near a Tricritical Point 

In the following we will discuss some relevant facts about 
the scaling near a tricritical point, see, e.g., Refl46l Consider 
some quantity O, that, in the thermodynamic limit, scales as 
0(£=~) ~ where £ is the correlation length. The Finite 
Size Scaling (FSS) ansatz, tells us how the same quantity be- 
haves in a finite system of size L. Close to the tricritical point 



at (p u T t = T c (p t )) 

0(L,p t + Sp,T t + ST) = L X G(W u T ,L y Pu p ) , (16) 

where G is a scaling function, and we have neglected scal- 
ing corrections. As stated in Eqs. Q and ([3]), there are two 
relevant scaling fields, the thermal field uj and the disorder 
field Up. Both uj and u p are functions of 8p and 8T, the de- 
viations from the tricritical point. If we work at uj = 0, we 
should expect that, at linear order, Kp|„ r _ x 8p. Then the 
phase transition is of the second order if 8p < 0, and of the 
first order if 8p > 0. 

Our main assumption will be that the Maxwell construc- 
tion, see Ref. 04] and the previous subsection, enforces the 
constrain u j — to an accuracy of order &(L~ D ) (this ex- 
pectation is well founded in the first-order part of the critical 
lineal). Hence, Eq. ( fTB") simplifies to 

0(L,p, Maxwell) = L x G(L y p(p-p t ))(l + 0(L yT - D )) . 

(17) 

So, the Maxwell construction allows us to employ standard 
FSS EH with an effective scaling-corrections exponent co = 
D — yj. The combination of Eqs. Q and ((3), standard RFIM 
scaling relations 15 and the numerical estimates in Ref. 1241 
yield m = d+ /3Rfim/ v Rfim _ L48 ( 2 ). 

A further irrelevant scaling field uq = 1 / log Q with expo- 
nent — 6 is also present. 13 Numerically, 9 = 1.468(2) 24 , while 
we expect CO = 1.48(2). These two exponents are so similar 
that, given our limited numerical accuracy, we shall not at- 
tempt to distinguish them. However, we remark that one ex- 
pects a larger amplitude of the scaling corrections for Q = 4, 
what is confirmed by our data (see Fig. [6J. 

D. Numerical Simulations and Thermalization Checks 

We considered concentration values 0.65 < p < 1 and lat- 
tice sizes 12 < L < 64. The precise values are indicated 
in Table [j] The p-resolution becomes denser close to the L- 
dependent position of the tricritical point. For all pairs (p, 
L) we simulated 500 samples, with the obvious exception of 
p = l. 

Each sample was simulated on a e-grid fine enough to allow 
for a correct spline interpolation, see Fig. [4] The simulations 
at the different e values were mutually independent. Hence, 
we faced an embarrassingly parallel computational problem, 
suitable for Ibercivis (with a caveat, see below). 

All samples were simulated for the same number of Monte 
Carlo steps, at every e value. However, the number of Monte 
Carlo steps did depend on e, as we explain now. First, we ran 
all samples at a given e-value for a fixed amount of Swendsen- 
Wang steps (e.g. 3 x 10 5 fori =64, or 2 x 10 5 fori =48), then 
we assessed thermalization. 

The thermalization check was the standard logarithmic 
data-binning: for any given value of e, we computed different 
estimates of the sample-averaged /3(e), using disjoint pieces 
of the Monte Carlo history. On the first bin, we included only 
the second-half of the Monte Carlo history (i.e., our safest data 
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FIG. 4. (Color online) In order to ascertain thermalization, we use 
the standard logarithmic data binning (data corresponding to j6(e), 
as computed for L = 64, p = 0.95). Bin 1 was computed from the 
sample-average of the last half of the Monte Carlo history on each 
sample (bin 2 corresponds to the second quarter of the Monte Carlo 
history, bin 3 to the second eighth, and so forth). Statistical compati- 
bility among the different bins is a strong thermalization check. Lines 
are cubic spline-interpolations for each bin. In order to demonstrate 
the importance of having a dense enough simulation grid (in partic- 
ular, close to high curvature regions), the spline interpolation in the 
blue line ignores the data at e = —0.2875. 



L 


Simulated p values 


12 


0.65, 0.675, 0.7, 0.725, 0.75, 0.775, 0.8, 0.825, 0.832, 0.85, 
0.875, 0.9, 0.925, 0.9375, 0.95 


16 


0.65, 0.675, 0.7, 0.725, 0.75, 0.775, 0.8, 0.825, 0.85, 0.854, 
0.875, 0.9, 0.925, 0.9375, 0.95 


24 


0.7, 0.725, 0.75, 0.775, 0.8, 0.825, 0.832, 0.845, 0.85, 0.875, 
0.9, 0.925. 


32 


0.75, 0.775, 0.8, 0.825, 0.85, 0.854, 0.8625, 0.875, 0.886, 
0.8875, 0.9, 0.925, 0.9375, 0.95, 0.975. 


48 


0.75, 0.775, 0.8, 0.825, 0.85, 0.8625, 0.875, 0.877, 0.8875, 0.9, 
0.925, 0.9375, 0.95. 


64 


0.8, 0.825, 0.85, 0.86875, 0.875, 0.8875, 0.9, 0.925, 0.9375, 
0.95. 



TABLE I. For each of the lattice sizes L, we indicate the values of p 
(the concentration of magnetic sites) for which we carried out sim- 
ulations. We shall need to regard the various quantities defined, as 
continuous functions of the density of magnetic sites, p. We shall 
need as well the corresponding p-derivatives. As a rule, we have ob- 
tained these functions of p through a cubic-spline interpolation of the 
data computed at these p-values. In fact, some of them were chosen 
in order to minimize the interpolation errors at some particularly im- 
portant values of p, see Tables [TT| and |III| Derivatives with respect to 
p were computed simply by derivating the cubic-spline interpolating 
function. The error estimates where obtained through a jack-knife 
(see for instance Ref. 46) over the sample-averages. 



from the point of view of thermalization). The second bin con- 
tained only the second quarter of the Monte Carlo history, etc. 
We checked for statistical compatibility, at least, among the 
first and second bins, see Fig. [4] If for a given value of e the 
thermalization criterion was not met, the total simulation time 
was doubled. The procedure was cycled until convergence 
was achieved. We note that, for the concentrations nearest to 
p = 1, we encountered strong metastabilities, that prevented 
us from simulating L— 128 (that could instead be simulated 
forg = 4inRef.[[l]>. 

The thermalization protocol is not well suited for Ibercivis, 
because the simulation of a given sample at some difficult en- 
ergy may last up to some days. Yet, Ibercivis relies on vol- 
unteers' computers that frequently switch from on-line to off- 
line. To minimize the number of unfinished simulations, we 
have implemented a continuity system. It divides every simu- 
lation, no matter how long it is, in small time steps (typically 
30 minutes). After every step, consistency checks are per- 
formed and the current system configuration is sent again to 
the simulation queue. This solved the problem for relatively 
long (5-6 hours) simulations but the few more demanding sim- 
ulations were completed on local clusters. 

Altogether, this work has consumed (the equivalent of) 3 x 
10 6 hours of a single Intel Core2 duo at 2.5 GHz. 

We should also mention that we have performed some new, 
short simulations for Q — 4 at p — 0.95, complementary to 
those reported in Ref. [19| The simulated sizes were L = 24 
and L = 48 (128 samples each). Our goal was to improve the 
accuracy of the interpolations described below. 



IV. RESULTS 

To check the Cardy-Jacobsen conjecture we have per- 
formed numerical simulations for 2 = 8, hence further ap- 
proaching the large-Q limit where the mapping becomes 

exactfSl 

Consider the p and L evolution of the latent-heat and the 
surface tension in Fig. [5] If p<pi (i.e. if we are in the second- 
order piece of the critical line), both Ae and E vanish in the 
large-L limit (the two are positive for p > p t ). However, for 
small lattices, both Ae and E decrease gently upon decreas- 
ing p which suggests that dilution merely smoothed the first- 
order transition. However, the curve for Ae becomes sharper 
upon increasing L. Indeed the Potts-RFIM mapping 13 implies 
Aeoc (p-p t )P with p=p KFm ~ 0.0l| P W hich is barely dis- 
tinguishable from a discontinuous jump. Furthermore, the L- 
dependent position of the tricritical point p\ (for instance, the 
point of sharpest drop of Ae in Fig. [5] — top) grows quickly 
with L. On the view of the D = 2 no-go theorems,^ one could 
be afraid that p\ — > 1 for large L also in D = 3. We know 
that this is not the casep but it is clear that a careful scaling 
analysis is needed. 

Eq. ( fT7| ) tells us that L 2 ~ e E is scale-invariant, and thus al- 
lows to locate the tricritical point (because X£ = — D + 1 = 
9 - 2, 9 = 1.469(20]!% Indeed, in Fig. |]— bottom, we see 
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FIG. 5. (Color online) Latent heat Ae (Top) and surface tension E 
(Middle) as a function spin concentration, p, for each lattice size 
(lines are linear interpolations). Lines end at the smallest p that al- 
lowed to perform the Maxwell construction for each L. Bottom: 
Scaled surface tension using 9 = 1.469(20pl(the lines joining the 
data are cubic splines). The vertical gray line shows the infinite vol- 
ume extrapolation for p t - 



that the curves for system sizes L\ < L 2 cross at p t 



Lite. 



Lf- E(L 1 ,p^)=L, 



(18) 



(Pt 



L,,Lo 



Pi when L\ diverges). We recall that a similar 



method was used recently in a spin-glass context. 53 There are 
two main consequences of choosing a wrong estimate of ex- 
ponent 9 in Fig. [5] — bottom and Eq. ([18). First, in the limit 
of large lattice sizes, the height of the crossing point diverges 
(or goes to zero) if 9 is underestimated (overestimated). Sec- 
ond, the size corrections to the crossing points are larger for a 
wrong 9. Specifically, p t l — pt = &(L { y ''). The amplitude 
for these scaling corrections cancels only for the exact choice 
of0. 

The critical exponent for a quantity O is obtained from its 



quotients at p t 



Lite {54 55 



0(L 2 ) 



0(Li) 



L X ,L 2 



*0 



l+A c 
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^2 



jw 



(19) 



Above, we included only the leading scaling-corrections (Ao 
is an amplitude). We use Eq. ( fT9) for the logarithmic p- 
derivative of E (scaling dimension x — y p ), and for the 
latent heat (scaling dimension x = f5y p , which should be 



JRF1M /..RFIM 



/v , according to Cardy and Jacobsen ). Our re- 
sults are in Table ^(Q = 8), and Table [TO] (g = 4). In both 



(12,16) 
(12,24) 
(16,24) 
(16,32) 
(24,32) 
(24,48) 
(32,48) 
(32,64) 
(48,64) 



Pt 

0.8947(38)(17) 

0.8942(16)(15) 

0.8939(28)(14) 

0.8966(13)(11) 

0.8989(28)(10) 

0.9031(14X10) 

0.9057(21X9) 

0.9040(H)(8) 

0.9026(21)(5) 



yp 

0.89(23)(2) 

0.82(g)(2) 

0.79(18)(3) 

0.85(13)(3) 

0.94(26)(3) 

0.80(6)(05) 

0.84(10)(1) 

0.86(5)(1) 

0.99(14)(3) 



0.108(5)(3) 
0.107(3)(3) 
0.107(6)(4) 
0.111(3X4) 
0.118(8X5) 
0.128(4)(6) 
0.139(8)(7) 
0.134(5)(7) 
0.126(10X8) 



py P 

0.095(9)(5) 
0.075(3)(4) 
0.061(5X3) 
0.050(2)(2) 
0.035(5)(2) 
0.027(2)(2) 
0.021(4)(2) 
0.023(3)(1) 
0.024(6)(1) 



TABLE II. Quotient-method for 2 = 8. For each pair of lattices 
(LijZ/j), we extract the crossing point pf'' Ll 



height of the crossing point, E(L\,pf ll ' La )I% 



see Eq. ( [18) , and the 
', The effective criti- 
cal exponents y p and f5y p are obtained using the quotients method, 
Eq. (19) . For each data, we indicate two error bars. The first 
error is statistical. The second error is due to the uncertainty in 
6 = 1.469(20)P 



(Li,L 2 ) 



Pt 



yp 



r 2 — 0^cross 



py P 



(16,24) 

(16,32) 

(24,32) 

(24,48) 

(32,48) 

(32,64) 

(48,64) 

(48,128) 

(64,128) 



0.9249(30X8) 

0.9324(19)(8) 

0.9400(30)(6) 

0.9455(19)(9) 

0.9506(27)(8) 

0.9489(13)(7) 

0.9473(31)(5) 

0.9491(9)(5) 

0.9497(14)(5) 



1.40(46)(3) 
1.11(20X5) 
1 .22(33)(1) 
0.83(8)(3) 

0.79(18)(3) 
0.78(9)(2) 

0.92(24)(3) 

0.77(8)(2) 

0.71(13)(2) 



0.0113(6X5) 

0.0125(5)(6) 

0.0159(12)(8) 

0.0179(9)(8) 

0.0215(17)(10) 

0.0206(9)(H) 

0.0191(25X12) 

0.0204(10X13) 

0.0213(17)(14) 



0.285(H)(6) 

0.230(6)(6) 

0.175(12)(4) 

0.135(5)(4) 

0.112(7)(3) 

0.095(4)(3) 

0.070(10)(3) 

0.048(4)(3) 

0.038(6X3) 



TABLE III. Quotient method for Q = 4 (data from Ref.[T9] improved 
through control variates and the addition of new runs near p t ). Same 
notations of TablelfTl 



cases we see that the convergence of p t 1 ' Ll to the thermody- 
namic limit is very fast. The height of the crossing point seems 
also stable with growing sizes. 

The results in Tables[ll]and[IIl]need to be extrapolated to the 
limit of infinite system sizes. This can be done by considering 
leading order scaling corrections, as in Eq. ( [19) . The extrapo- 
lation greatly improves by imposing to Q = 4 and 8 a common 
extrapolation and the same scaling-corrections exponent co, as 
required by the Universality predicted in Ref.[T3l 

In this way, we obtain j8v p = 0.0022(48) (3) and co = 
1.36(8) (1) where the second parenthesis indicates the uncer- 
tainty induced by the error in ePThe fit quality is assessed 
through the % 2 test. We obtain £ 2 /dof = 8.5/14, which is al- 
most too good (dof stands for the number of degrees of free- 
dom of the fit). Indeed the probability of getting such a low 
value of x 2 with 14 degrees of freedom is only 14%. We note 
as well that j5y p = 0.0022(48) (3) is only barely compatible 
with the best RFIM estimate f5y p = 0.01 19(4)53 (since the 
discrepancy is as large as two standard deviations). 
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0.01 0.02 
(L" 1 "-L 2 w )/log(L 2 /L 1 ) 

FIG. 6. (Color online) Determination of the scaling correction ex- 
ponent ft), from the size-dependent effective exponent fiy p for the 
vanishing latent-heat at the tricritical point as computed with the 
quotients method, see Tables [II] and [III] A common extrapolation 
Py p = 0.01 19(4 jSlis imposed in the joint-fit for the g = 4 and 2 = 8 
data. The figure of merit x 2 
variance matrix. 



14/15 is computed with the full co- 



At this point, we can try to disproof universality. We make 
the assumption that j5y p takes exactly the RFIM value, redo 
the fit and see the outcome of the % 2 test. This second fit, 
with Py p fixed to 0.01 19, turns out to be perfectly reasonable 
(£ 2 /dof = 14/15, see Fig. |6|. Hence we conclude that our 
data set is statistically compatible with universality. 

A second, unexpected bonus of fixing /5y p in the fit to 
the RFIM value, is a remarkable increase in the accuracy of 
CO = 1.53(5)(3), in excellent agreement with our expected 
CO = 1 .4 8(2) (r emember that co = D-y T = 9 + ^ RFlM /v mM , 
see Sec. |IIIC| >. Furthermore, from this value of co, we obtain 
Vt = D co = 1.47(8), in nice agreement with the large-g 
computation V7 = 1.49(9)p2l 

Following the same approach for y p , which is expected to 
coincide with l/v* 1 ™, we obtain y p = 0.775(46)(1) if we 
impose CO = 1.36 (we get y p = 0.779(41)(1) by taking CO = 
1.53). Both fits are fair (% 2 /dof = 13.6/15 and X 2 /dof = 
13.2/15, respectively). 

Our y p is in the lower range of previous numerical and an- 
alytical estimates: 0.73 < l/v RF1M < 1.12p3E3SlHyperscal- 
ing and our y p implies a slightly positive specific heat expo- 
nent a = (2y p —D + 9)/y p — 0.03(10), in agreement with 
experimental claims of a (possibly logarithmic) divergence.^ 
We warn however that severe hyperscaling violations [namely 
a = —0.63(7)] have been reported in numerical work!^ 

One may compute as well the exponent 9, by fitting 
E (L, p[' 2L ) =A Q L e (l+ B Q L- m ) (only the amplitudes A Q and 
Bq are Q-dependent on the fit). Taking CO = 1.5(1), we ob- 
tain = 1.52(H)(2), with an acceptable fit (x 2 /dof = 4.9 /3). 
The result is compatible with, but less accurate than, the latest 
RFIM result 9 = 1.469(20)P 



V. CONCLUSIONS 



In summary, we have presented a finite-size scaling analysis 
of the tricritical point of the site-diluted Potts model in three 
dimensions for 2 = 4 and 8 internal states. By considering 
leading-order scaling corrections we were able to show that 
the relevant Universality class for the tricritical point is the one 
of the RFIM. To our knowledge, this is the first verification of 
the Cardy-Jacobsen conjecture!^ 

Three technical ingredients were crucial to obtain this 
achievement: the use of the microcanonical Monte CarloJSl 
a new definition of the disorder average, 19 and the use of the 
citizen supercomputer Ibercivis!^ 
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Appendix A: Control variates 

The statistical quality of data may sometimes be signifi- 
cantly increased by means of a very simple trick, named con- 
trol variates (see e.g. Ref.l58l. 

In short, we want to improve our estimation of a stochastic 
variable A through its correlations with another random vari- 
able B (B is named a control variate). If B — and A = A + aB, 
then the expectation value does not change: A = A. How- 



0.860686 



0.860684 



0.860682 




0.949 



FIG. 7. (Color online) Scatter plot of each sample's inverse critical 
temperature vs. the concentration of magnetic sites, £,•£,•/¥. Data 
for 500 samples of L = 64 and p = 0.95. The correlation coefficient 
that gives the optimal coupling to the control variate, see Eq. |ATJ, is 
a* = 0.956. 
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FIG. 8. (Color online) Inverse critical temperature j3 C £ as a func- 
tion of an inverse lattice volume, 1/L 3 . Data obtained for p = 0.95. 
The error reduction obtained with control variates is significant (blue 
points). In the linear fit we considered only data with L > 16 and 
improved through control variates. 
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FIG. 9. (Color online) The size dependent latent-heat for p = 0.95, 
as a function of the inverse lattice volume. The linear fit includes 
only data with L > 16 that were improved through control variates. 



ever, depending on the arbitrary election of a, we can get 
var(A) < var(A). The a election minimizing the variance 
var(A) is 



^/var(A)var(Z?) ' 



which coincides with the correlation coefficient r^g. The op- 
timal variance is 



var(A*)=var(A)(l-ri fi ). (A2) 



Hence, the stronger the statistical correlation (or anticorrela- 
tion) between A and B, the more effective the control variate 
is. 

In our case, a rather obvious control variate is 



namely the difference among the real and the nominal concen- 
trations of magnetic sites. It is clear that the disorder average 
B vanishes. We will employ B to improve the determination 
of the sample-averaged /3(e). Note that, although the value of 
B does not depend on the considered energy (it is fixed by the 
{£,■}), its correlation coefficient with (j5) e needs to be com- 
puted for all energies in the e-grid. 

B is extremely effective as a control variate for the compu- 
tation of the inverse critical temperature j3 c , as suggested from 
Fig. [7] The correlation coefficient in that plot is so high, 0.956, 
that the expected error reduction factor is 3.4. However, the 
alert reader will note that this is a hasty conclusion. In fact, 
the j3 c obtained from /3(e) is not exactly the average of the in- 
verse critical temperatures found for each sample. The reason 
for this non-linearity in the Maxwell rule, see Eq. ( [TT| ), is that 
the energies e A o are not the same for j3 (e) and for the (j3) e in 
a given sample. Yet, the dependency on e A ° of the integral in 
Eq. ( fTTj ) is extremely weak [recall the stationarity condition 
with respect to e in Eq. p0)]. 

In fact, the correct computation with /3 (e) does show a sig- 
nificant error reduction, see Fig. [8] close to the factor 3.4 an- 
ticipated by the naive analysis in Fig. [7] We note in Fig. [9] 
an equally significant reduction of the statistical errors for the 
latent heat. Therefore, our computation of these quantities, 
obtained with only 500 samples, has been made equivalent to 
a 5000-samples computation. This is a remarkable reward for 
such a simple analysis. 
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